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An acoustic analysis based on the Ffowcs Williams 
and Hawkings equation was performed for a high-lift 
system. As input, the acoustic analysis used un- 
steady flow data obtained from a highly resolved, 
time-dependent, Reynolds- averaged Navier-Stokes 
calculation. The analysis strongly suggests that vor- 
tex shedding from the trailing edge of the slat results 
in a high-amplitude, high-frequency acoustic signal, 
similar to that which was observed in a correspond- 
ing experimental study of the high-lift system. 

Introduction 

Airframe-generated noise is an important compo- 
nent of the total noise radiated from commercial 
aircraft, especially during the approach portion of 
the flight path when the engines are run at reduced 
power. Recent studies by Davy and Remy 1 on a 
scale model Airbus aircraft indicate that the high- 
lift devices and landing gear are the main sources 
of airframe noise when the aircraft is configured 
for approach. Earlier tests on a model of a DC- 
10 also identified the high-lift system as an impor- 
tant airframe noise source. 2 Dobrzynski, Nagakura, 
Gehlhar, and Buschbaum 3 performed full-scale ex- 
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perimental studies in an open-jet wind tunnel of a 
portion of a wing equipped with a high-lift system. 
They found that both the leading-edge slat and the 
side edge of the trailing flap contributed significantly 
to the airframe noise. 

An extensive experimental and computational ef- 
fort to study the various mechanisms associated 
with airframe-generated noise is currently underway 
at NASA Langley Research Center. 4 Considerable 
progress has already been made in understanding 
various aspects of the noise-generation process on 
the flap side edge. 5-11 More recently, attention has 
focussed on noise generated in the vicinity of the 
leading-edge slat. 

A cooperative test involving NASA’s High-Lift 
Program Element and NASA’s Airframe Noise Team 
was conducted in NASA Langley Research Center’s 
Low- Turbulence Pressure Tunnel (LTPT). Variation 
of the pressure in the tunnel allows the Reynolds 
number to be changed at constant Mach number. 
In these tests, the Reynolds number based on the 
cruise-wing chord varied from 3.6 to 19 million. No 
qualitative changes were observed in the data for 
Reynolds numbers above 7.2 million. The model 
tested in the tunnel is known as the Energy Effi- 
cient Transport (EET) model. 12 The EET model 
tested includes a full-span leading-edge slat and a 
part-span trailing flap. To obtain acoustic data, 
members of Boeing Commercial Airplane Company 
designed and built a microphone array that was in- 
stalled in the ceiling of the wind tunnel. The mi- 
crophone array and the subsequent data processing 
followed techniques developed earlier at Boeing. 13 
These techniques have previously been used suc- 
cessfully to determine the noise radiated from local- 
ized sources, even in hard-walled, non-anechoic wind 
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Figure 1. Schematic of model in wind tunnel. Flow 
from left to right. Acoustic directivity angles 0 and 
270 deg. are indicated. View is rotated relative to 
experimental setup. 


Figure 2. Acoustic spectrum based upon 1 / 12th oc- 
tave bins with array focussed on slat region. Config- 
uration angle of attack is 10 deg., Reynolds number 
is 7.2 million, Mach number is 0.2. 


tunnels like the LTPT. To measure noise radiating 
groundward in normal flight, the EET model was 
mounted upside-down in the tunnel so that the pres- 
sure surface faced the array in the wind-tunnel ceil- 
ing. Flowever, to reduce confusion, all references to 
directions in this paper will conform to the schematic 
of the experimental setup shown in Fig. 1. The view 
in Fig. 1 has been rotated from the physical orien- 
tation to one that is more intuitive. An enormous 
amount of data was collected and continues to be 
analyzed. 

Figure 2 illustrates one unexpected and until re- 
cently, perplexing feature of the experimental data. 
For the case in which the slat deflection, 6 S , is 30 
degrees, a very large amplitude peak is observed in 
the acoustic spectrum in the vicinity of 50 kFlz. This 
peak rises almost 20 dB above the signal observed 
for the case in which the slat is deflected 20 degrees. 
During the course of the experiment, efforts to elim- 
inate the high-frequency peak by altering the over- 
hang of the slat were largely unsuccessful. Only for 
cases in which the overhang became unrealistically 
large was a significant change in the high-frequency 
acoustic peak observed. Increasing the configura- 
tion’s angle-of-attack from 10 to 15 degrees, reduced 
the amplitude of the high-frequency peak by approx- 
imately 10 dB. For some time, no consistent explana- 
tion of the observed phenomena was available. The 
focus of this paper and a companion paper 14 is to ex- 
plain the observed large-amplitude, high-frequency 
peak in the experimentally obtained acoustic spec- 


trum. 

Khorrami et. al 14 provides details of unsteady, two- 
dimensional (2D), Reynolds-averaged Navier-Stokes 
(RANS) calculations designed to mimic the exper- 
imental conditions. Ri particular, the RANS com- 
putation was specially designed to properly incor- 
porate and resolve the small, but finite t. railing-edge 
thickness of the slat. Extremely small grid cells were 
used in the vicinity of the slat trailing edge and the 
time step was chosen to ensure more than 120 time 
steps per period of a 50 kHz signal. Initially calcula- 
tions were performed with a slat trailing edge thick- 
ness h of approximately 0.07 percent of the cruise- 
wing chordlength C , which was an estimate of the 
actual slat t. railing-edge thickness. Slat deflections 
of both 30 and 20 degr<$ls were simulated. These 
calculations clearly show vortex shedding from the 
slat trailing edge for the case with a 30 degree slat 
deflection. Figure 3 shows a snapshot of the pressure 
fluctuations produced in the flow held. For physical 
reasons that are not yet clearly understood, the vor- 
tex shedding virtually disappears for the case of a 
20 degree slat deflection. Acoustic analyses of the 
first set of data suggested that the initial grid dis- 
tribution was insufficient to completely resolve the 
complex acoustic held in the cove and other regions 
in the vicinity of the slat. Later calculations with 
an enhanced grid were performed for a 30 degree 
slat deflection with h/C 0.0009. The slat thick- 
ness was adjusted to more accurately represent the 
measured slat thickness on the model. In all cases, 
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Figure 3. Instantaneous fluctuation pressure, in 
vicinity of leading-edge slat, from CFD calculation. 
Slat deflection is 30 deg. Wiggles at edges of dark 
and light, bands are contouring artifacts. 

the high cost of the calculations limited the duration 
of the temporal sample that was obtained. Here we 
discuss the use of that limited sample of unsteady 
computational data to perform acoustic analyses of 
the generated noise. 

Acoustic Procedure 

Previously Singer et. a.l 15 explored the use of un- 
steady computational data, in acoustic-propagation 
codes based on the Ffowcs Williams and Hawk- 
ings 113 (hereafter referred to as FW-H) equation. 
Such codes compute the acoustic signal at. a. dis- 
tant. observer position by integrating the FW-H 
equation. Following Brent.ner and Fa.ra.ssa.t., 1 ' the 
FW-H equation may be written in differential form 
as 

° V(x ’ i)= a^7 [l;jW)1 

-j£iw(/)j+ |tws,xw)j a) 

where: n 2 = — V 2 is the wave operator, c is 

the ambient, speed of sound, t is observer time, p' is 
the acoustic pressure, p' is the perturbation density, 
po is the free-st.rea.m density, / = 0 describes the 
integration surface, 6(f) is the Dirac delta, function, 
and H(f) is the Heaviside function. The quantities 


Ui and Li are defined as 

Ui — (1 )'l’i H (2) 

Po Po 

and 

Li = I’ n + pui(u n - v n ) (3) 

respectively. In the above equations, p is the total 
density, ptq, is the momentum in the i direction, 
i>i is the velocity of the integration surface / = 0, 
and Pij is the compressive stress tensor. For an in- 
viscid fluid, Pij = p'Sij where &ij is the Kronecker 
delta. The subscript, n indicates the projection of a. 
vector quantity in the surface-normal direction. To 
obtain a. solution to Eq. (1), the first, term on the 
right-hand-side must, be integrated over the volume 
outside the integration surface / = 0 wherever the 
Light. hill stress tensor Tij is nonzero in this region. 
In the work reported here, this term is neglected. 
However the main effects of nonzero Tij within the 
flow can be included by choosing an integration sur- 
face that, contains all of the volume with significant. 
Tij contributions. 

The other terms on the right-hand-side of Eq. (1) 
include terms that, are determined by the unsteady 
flow-field data, on the integration surface. Provided 
that, the unsteady flow data, on the integration sur- 
face / = 0 is correct., Reference 15 demonstrated 
that, the FW-H equation correctly propagates the 
acoustic radiation from several source regions, in- 
cluding the complex signals associated with acoustic 
scattering from sharp edges. 

The extremely small time step required in the 
RANS calculation to adequately resolve the high- 
frequency flow physics resulted in a. limited tempo- 
ral duration of the data. The total time represented 
by the unsteady calculation was 0.68 ms, approxi- 
mately long enough for an acoustic wave generated 
at. the leading edge of the main element, to propa- 
gate halfway down the cruise- wing chord C . Because 
of the short, time duration, the initial transient, is 
important, even when attention is restricted to the 
vicinity of the slat. Figure 4 shows the perturba- 
tion pressure signal at. the slat, trailing edge, the slat, 
cusp, and the slat, leading edge. To account, for the 
rapidly decaying amplitude of the fluctuating pres- 
sure, relative to the scale of the pressure fluctuations 
at. the slat, trailing edge, the scale of the pressure 
fluctuations is magnified by a. factor of 40 at, the slat, 
cusp and a, factor of 200 at the slat, leading edge. 
Although the data, at, the slat, trailing edge is qua.si- 
periodic from the start,, over 0.2 ms pass before the 
transient, goes by the slat, cusp and almost, 0.38 ms 
pass before the transient, propagates past, the slat, 
leading edge. After passage of the initial transient,, 
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Figure 4. Perturbation pressure as function of time. 

Slat deflection is 30 degrees. slat trailing 

edge, slat cusp, slat leading 

edge. Vertical scale for data at slat cusp is mag- 
nified by factor of 40 and shifted. Vertical scale for 
data at slat leading edge is magnified by a factor of 
200 and shifted. 

a quasi-periodic condition prevails at all locations. 
To limit the effect of the transient, the first. 25 per- 
cent of the data record is not used to produce the 
results presented here. Auxiliary calculations sug- 
gest that the use or nonuse of the first. 25 percent, 
of the data, record produce relatively small quanti- 
tative variations in the results. All of the qualitative 
results are unchanged. To preclude any confusion, 
hereafter the term “data, record” or any equivalents 
should be construed to refer to only the portion of 
the data. that, is actually used. 

To deal with the limited data, time series, a. mod- 
ified Planning window is applied to the data. The 
modified Planning window includes a. standard Plan- 
ning filter for the first, and last. 12.5 percent, of the 
data, and a. boxcar filter for the middle 75 percent.. 
The windowed data, is scaled to preserve the original 
energy in the signal. The resulting data, sequence is 
then implicitly repeated as needed to provide an in- 
put. signal of arbitrary duration. The windowing cre- 
ates an artificial periodicity at. approximately 1960 
Hz, but. because this frequency is much lower than 
the vortex shedding frequency, the artificial period- 
icity does not. introduce any problems. 

The application of acoustic theories to 2D flow 
duly is a. problem that, is likely to become more visi- 
ble as computational fluid dynamics (CFD) is relied 
upon more regularly to provide unsteady flow data. 


for use in acoustic calculations. Time-dependent., 
3D, CFD data, is extremely expensive to produce. 
In many applications, as in the current, problem, 
the primary aerodynamic phenomena, that, generate 
noise are essentially 2D. The 3D effects are largely 
restricted to the fact. that, the true 3D unsteady flow 
structure is not. completely correlated in the third 
direction. Cox et. at. 18 computed 2D and 3D vortex 
shedding over circular cylinders and then used the 
results to calculate the acoustic pressure at. an ob- 
server location. To use the 2D CFD computations, 
the acoustic calculations were performed by assum- 
ing perfect, spa.nwise correlation of the flow over a. 
finite span. They found that, the acoustic ampli- 
tude could increase by as much as 20 dB simply by 
extending the finite span from 5 to 100 cylinder di- 
ameters. 

As a. first approximation, a. 2D version of the FW- 
H equations is used to predict, the sound held. Here 
we use the code developed by Lockard 19 for com- 
puting the 2D acoustic held from 2D CFD data. 
As noted above, we expect, the 2D results to have 
greater amplitudes than those observed in the exper- 
iment., but. the qualitative features of the acoustics 
are not. expected to differ substantially. To study the 
effects of the spa.nwise correlation length, a. limited 
number of 3D acoustic calculations are performed 
using the same FW-H code as was used in Ref. 15. 
As input, to the 3D code, the 2D CFD data, is re- 
peated for a. finite distance in the spa.nwise direc- 
tion. Previous tests on idealized problems confirm 
that the 3D FW-H code and the 2D FW-H code 
give identical results for model problems when the 
spa.nwise extent, is sufficiently long. 

For consistency, all of the acoustic calculations are 
performed for observers located a. hxed distance from 
the trailing edge of the slat. The hxed distance cor- 
responds to the distance from the slat, trailing edge 
to the centroid of the acoustic array. Directivity 
angles are indicated in Fig. 1; 0 degrees is in the 
downstream direction, 270 degrees is groundward in 
normal flight., towards the microphone array in the 
wind tunnel. 

Another important, issue involves the choice of in- 
tegration surface for the FW-H calculation. Figure 5 
illustrates the two integration surfaces that, have 
been used for the FW-H calculations for the cases 
with a. slat t.ra.iling-edge thickness h/C Pri 0.0007. 
The solid lines correspond to the component, surfaces 
that, lie on the solid bodies of the slat, and main ele- 
ment.. This combination of surfaces is designated the 
“on-body surface.” Because the limited time sam- 
ple is insufficient, for acoustic signals to propagate 
from the leading-edge slat, to the- trailing flap, the 
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Figure 5. Integration surfaces used for FW-H cal- 
culations with h/C & 0.0007; on-body sur- 
face; off-body surface. 

information on the flap is not included. If it wer# 
included, it would contribute nothing to the noise 
because the computations show that the flow in the 
vicinity of the flap is steady. The dotted line cor- 
responds to a permeable mathematical surface that 
contains within it the boundary layers on all t.hreq 
elements of the high-lift, system as well as the re- 
gion where the shed vortices reside. This integration 
surface includes all the major noise sources in this 
flow, and is designated the “off-body surface.” Be- 
cause the flow-field quadrupole contributions to the 
noise are contained within the off-body surface, i lie 
neglected term in Eq. (1) is essentially zero for 
this case and little approximation is made by not in- 
tegrating it explicitly. Although thig; consideration 
suggests that the off-body surface is a better choice 
than the on-body surface, in practice CFD grid res- 
olution issues complicate the choice. Because of the 
limited order of accuracy of the CFD calculations, 14 
flow fluctuations that are present on the on-body 
surface may have been numerically dissipated and 
dispersed by the CFD computation before the fluc- 
tuations arrive on the off-body surface, especially in 
the region between the pressure surface of the slat 
and the leading edge of the main element. 

The deficiencies in the CFD grid between the slat 
and the main element have been minimized by refin- 
ing the grid for the calculation with the slat trailing- 
edge thickness h/C Pri 0.0009. In addition, our ex- 
perience with the prior calculations suggested that 
data be saved on a number of additional integration 
surfaces. Close-up views of these integration sur- 
faces in the vicinity of the slat are shown in Figs. 6 
and 7. For convenience, the four integration sur- 
faces used for the h/C Pd 0.0009 case are labeled 
0-3. A brief description of each surface is given in 
Table 1. In Fig. 6 the solid lines (surface 0) repre- 
sent the physical surfaces and correspond to the on- 
body integration surfaces used in the h/C Pd 0.0007 
cases. The dotted lines (surface 3) indicate a surface 
that includes all boundary layers and wakes and cor- 



Figure 6. Integration surfaces used for FW-H cal- 
culations with h/C Pd 0.0009; on-body (sur- 
face 0); outside of boundary layers and 

cove (surfaces 1 and 2, see Fig. 7 to distinguish); 

far outside of wakes and boundary layers 

(surface 3). 

responds closely to the off-body integration surface 
used in the h/C Pd 0.0007 cases. The dashed lines 
(surfaces 1 and 2) are intermediate surfaces that ex- 
tend just outside of the boundary layers and include 
the entire cove region between the slat and main el- 
ements. The distinction between surface 1 and 2 is 
shown in Fig. 7. Surface 1 (shown with the dashed 
line) does not include the shed vortices in the wake of 
the slat while surface 2 (shown with the dotted line) 
includes a region that contains a limited number of 
the shed vortices. Surfaces 1 and 2 are believed to 
offer the most reasonable compromise between in- 
clusion of acoustic sources and adequate CFD grid 
resolution from the sources to the integration sur- 
faces. 

Results 

Figure 8 shows computed spectra based on 1 /12 th 
octave bins for an observer located at 270 degrees. 
Only results for the on-body integration surfaces are 
shown. For the slat deflection of 20 degrees, the 
small-amplitude peak around 45 kHz is a residue of 
the initial transient. The analysis of the data in 
Ref. 14 shows no sustained vortex shedding in i lie 
CFD flow results for the 20 degree slat deflection. 
The corresponding FW-H spectrum essentially con- 
firms the intuitively obvious result that no signifi- 
cant noise radiates. On the other hand, the vortex 
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Figure 7. Close-up view of slat trailing edge;-. 

surface 0; surface 1; 

surface 2. Note that surface 1 and surface 2 share 
boundaries away from slat trailing edge. 


Surface Desciption 
n I on body 


off body - just, outside 
boundary layers and cove 
includes surface 1 and 
small region behind 
slat trailing edge 
far outside of all wakes 
and boundary layers 


Table 1. Description of integration surfaces. 



frequency, kHz 

Figure 8. Spectra for observer positioned at 270 deg 
using on-body integration surfaces; 30 deg 


slat deflection, h/C PS 0.0007; 20 deg slat 

deflection, h/C Pd 0.0007; 30 deg slat 

deflection, h = 0.0197 in. 


shedding associated with the 30 degree slat deflec- 
tion produces intense noise, peaking around 45 kFlz 
for the h/C Pd 0.0007 case and at about. 38 kFlz 
for the h/C Pd 0.0009 case. As expected, the thicker 
slat, trailing edge results in stronger disturbances and 
louder noise levels. The lower frequency peak for the 
thicker trailing edge is consistent, with the expec- 
tation of St.rouha.l scaling for the vortex shedding. 
Flowever, the 27%) increase in t.ra.iling-edge thickness 
results in an approximately 16%) decrease in domi- 
nant. frequency. Clearly, the t.ra.iling-edge thickness 
is not. the only important, flow parameter. Hammond 
and Redekopp 20 suggest, the use of the body thick- 
ness plus the displacement, thicknesses as the appro- 
priate length scale for St.rouha.l scaling. Figure 9 
shows the spectra, for both t.ra.iling-edge thicknesses 
when comparable off-body integration surfaces are 
used. The dominant, frequencies and general treads 
remain unchanged. 

Although the frequency of maximum noise does 
not. vary with differing integration surfaces, the di- 
rectivity pattern is sensitive to the choice of integra- 
tion surface. Figure 10 shows directivity patterns 
for the on- and off-body integration surfaces with 
h/C PS 0.0007. In this figure, at. any angle, the 
distance from the origin is proportional to the inte- 
grated pressure spectrum that, radiates at the chosen 
angle. The solid line represents the directivity com- 
puted with CFD data, on the on-body surface while 
the dashed line represents the directivity computed 
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Figure 9. Spectra for observer positioned at 270 
deg using off-body integration surfaces with 30 deg 

slat deflection; h/C Pd 0.0007; 

h/C Pd 0.0009. 

with the CFD data from the off-body surface. Signif- 
icant qualitative differences exist, especially with re- 
gard to the strong upstream and downward directiv- 
ity projections associated with the on-body integra- 
tion surface. Without experimental directivity mea- 
surements or another calculation that better resolves 
the acoustic signals in the region between the on- 
and off-body integration surfaces, choosing which 
integration surface provides the more realistic pre- 
diction is difficult. Fortunately, the h/C Pd 0.0009 
calculation provides adequate grid resolution to a 
variety of off-body surfaces. 

Figure 11 shows the directivity patterns that are 
obtained from the integration surfaces used with 
the thicker trailing-edge case. Many features of 
the directivity patterns are similiar for both the 
h/C Pd 0.0007 calculation and the h/C Pd 0.0009 cal- 
culation. In particular, for the on-body integration 
surfaces, strong noise radiation is directed slightly 
above and below 180 degrees and at approximately 
90, 225, and 270-300 degrees. These strong directiv- 
ity projections exist for both slat trailing-edge thick- 
nesses. As would be expected, the presense of the 
main element blocks significant acoustic radiation in 
the downstream direction. The use of off-body inte- 
gration surfaces has a similar effect for both trailing- 
edge thickness cases. The strong upstream-directed 
projections slightly above and below 180 degrees, 
which are associated with the on-body integration 
surfaces, merge and have reduced levels. The acous- 
tic radiation peak directed at approximately 225 de- 



p'cos(e) 

Figure 10. Directivity of acoustic signal for h/C Pd 

0. 0007; on-body surface; off- 

body surface. 

grees essentially disappears when off-body integra- 
tion surfaces are used and the strong downwardly- 
directed acoustic radiation is reduced in amplitude. 

The directivity patterns for off-body surfaces 1- 
3, shown in Fig. 11 are generally similar. How- 
ever, where differences occur, simple explanations 
or trends are lacking. For instancfte surface 3, which 
includes the entire slat wake, suggests stronger ra- 
diation in the 90-120 degree range than does either 
surface 1 or 2. However, the use of surface 2, which 
includes some of the slat wake, results in less acoustic 
radiation in this region than does the use of surface 

1, which includes none of the slat wake. Therefore, 
no clear trend exists with respect to the inclusion of 
additional portions of the slat wake. Similarly, no 
clear trend exists in the 330-360 degree range. In 
this angular range, the noise radiated from surface 
3 is less than that radiated from either surface 1 or 
or surface 2, but the noise radiated from surface 2 is 
slightly greater than that radiated from surface 1. 

Although the acoustic array used in the exper- 
iment was not intended to provide any directiv- 
ity information, the high-frequency acoustic signal 
was so loud that it overwhelmed the intrinsic wind- 
tunnel noise and can be identified from the spec- 
trum of some of the individual microphones used 
in the acoustic array. The relative amplitudes of 
the mean square fluctuating pressure in a frequency 
range around 50 kHz from a subset of microphones 
having approximately the same cross-stream loca- 
tion are shown in Fig. 12. The abscissa in the figure 
indicates streamwise distance. The microphone lo- 
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Figure 11. Directivity of acoustic signal for 

h/C & 0.0009; on-body surface (surface 0); 

surface 1 ; — surface 2 ; 

surface 3. 

cations are shown as squares in the figure and their 
positions relative to the airfoil are easily deduced. 
The computed mean square fluctuating pressure pro- 
jected to the acoustic-array surface is shown as the 
dashed line. The computed results are calculated 
using integration surface 1. The maximum ampli- 
tude of the microphone data is scaled with the max- 
imum amplitude of the calculation. Far upstream of 
the airfoil, both the microphone response and the 
computed noise level are flat. The non-zero mi- 
crophone response is probably associated with the 
wall-pressure fluctuations of the turbulent boundary 
layer along the wind-tunnel ceiling. These fluctua- 
tions are not included as part of the CFD calcula- 
tions. Slightly upstream of the slat leading edge, the 
noise level rises. Unfortunately the construction of 
the wind tunnel prevented the positioning of micro- 
phones over an extensive region that would include 
the streamwise location of the slat trailing edge. The 
maximum amplitude occurs in the mid-chord region 
followed by a sharp drop in amplitude. The qualita- 
tive features of the computations agree remarkably 
well with the microphone data and have been used to 
aid in the redesign of the acoustic array for a future 
wind-tunnel tipt. 

In the actual flow, the vortex-shedding process is 
not perfectly correlated in the spanwise direction. 
The correlation between two points is reduced as 
the spanwise distance between them increases. The 
physical decorrelation process could be modeled as a 
random walk, however numerical integration over a 


Figure 12. Comparison of squared acoustic pressure 
at individual microphones to that predicted com- 
putationally for h/C & 0.0009 case using integra- 
tion surface 1. Microphone positions and values are 
shown with squares; dashed line indicates computa- 
tionally predicted values. 

very large distance in the spanwise direction would 
then be required. Here we model the effects of the 
span by assuming perfect correlation along a fixed 
spanwise distance with no contribution from points 
outside of this distance. The half-span used in the 
calculations is denoted by L z and it is normalized 
with the cruise- wing chord C. The observers are 
positioned a distance d, = 2C from the slat trailing 
edge. 

To estimate the sort of behavior that we might ex- 
pect for the 3D acoustic calculation, we first exam- 
ine the analytically predicted behavior of a contin- 
uous line source in the spanwise direction. For this 
source, the contribution to the pressure will vary in 
the spanwise direction as 

h l . . exp (iwd/c^Jl + (z/d) 2 ) (4) 

where z is the spanwise coordinate and the ob- 
server is located at z = 0. The leading factor 
1/^/1 + (z/d.) 2 controls the decay of the amplitude 
of the acoustic contribution as the source points 
move away from the observer. The terms in the 
exponential determine the spatial periodicity. For 
small z, the spanwise oscillation period is large, but 
the period decreases with increasing z. To reduce 
the spanwise oscillation amplitude to 10 percent of 
its z = 0 value requires that z Pd 10 d, or, for the 
cases studied here, about. 11 meters. This distance 
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Figure 13. Directivity of acoustic signals for h/C Pri 

0.0009; L.JC = 0.25; L z /C = 

2.0; 2D calculation. 

is far greater than the wind-tunnel span, and would 
clearly be much greater than any reasonable esti- 
mate of the spanwise coherence. These estimates 
suggest that the 2D calculations are likely to sub- 
stantially overpredict. the amplitude of the acoustic 
field. Further analysis of Eq. (4) suggests that near 
z = 0, the spanwise oscillation period z p for a 38 kFlz 
source is z p /C Prf 0.26, but that for large z, the os- 
cillation period is much smaller, with z p /C Prf 0.015. 
The variation of the oscillation period requires that 
spanwise grid-resolution requirements be reviewed 
whenever the spanwise extent changes. 

Because 3D acoustic calculations are much more 
expensive than their 2D counterparts, 36 observers 
(spaced ten degrees apart.) are used in an initial set 
of 3D acoustic calculations to determine trends with 
variation of the spanwise extent L z . Only data on in- 
tegration surface 1 was used for these calculations. 
After ascertaining the spanwise resolution require- 
ments and finding no unexpected tendencies, acous- 
tic calculations using L z /C = 0.25 and L z /C = 2.0 
were performed with 180 observers spaced two de- 
grees apart.. The directivity patterns are compared 
with the results from the 2D calculation in Fig. 13. 

Examination of Fig. 13 reveals some interesting 
trends. As expected, with increased spanwise extent, 
the noise increases. Cases run with 0.25 < L z /C < 
2.0 show that the increase in noise is not. monotonic, 
varying instead with the spanwise period. Cox et. 
a.l. 18 observed a. similar trend until, at some value 
of L z , the average noise level a.sympt.ot.ed and only 
small variations associated with the periodicity in z 



Figure 14. Directivity of acoustic signals for h/C Pri 

0.0009; L z /C = 2.0 scaled by a. factor of 1.4; 

2D calculation. 


were observed. Our data, suggest, that the qualita- 
tive features of even a. complex directivity pattern 
establish themselves at relatively small values of L z . 
In Fig. 13 all of the directivity lobes found in the 2D 
results can be traced back to corresponding lobes 
obtained with both L z /C = 2.0 and L z /C = 0.25. 
Recall that L z /C Prf 0.26 corresponds to the span- 
wise period in Eq. (4) for z Pd 0; therefore the basic 
directivity features are established within the first, 
spanwise period. 

Although uniform scaling of the sma.ll-t.o- 
modera.t.e L z results does not. exactly replicate the 
2D directivity pattern, such a. simple scaling is a. 
reasonable estimate. A comparison of the 2D re- 
sults and the L z /C = 2.0 results scaled by a. factor 
of 1.4 is shown in Fig. 14. The agreement, between 
the two is surprisingly good. However, the scaling 
factor was determined empirically. Efforts to derive 
the scaling factor have been unsuccessful thus far. 
For instance, the ratio obtained by dividing the am- 
plitude of the integration of Eq. (4) from z = — oo 
to z = oo by the amplitude of the: the integration 
from z = — L z t.o z = L z might, be expected to ap- 
proximate the scaling factor. However, this ratio is 
in the range of 1.03, rather than the observed 1.4. 
Variations on this theme, such as taking ratios of 
only the real or imaginary parts produced similarly 
disappointing results. Additional work still needs to 
be done to take full advantage of the 2D acoustic 
solution in this type of qua.si-2D flow. 
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Conclusions 

Data from highly resolved, unsteady, 2D RANS 
calculations of a high-lift configuration with a 
blunted slat trailing edge were used to per- 
form acoustic analyses using the permeable-surface 
FW-H equation. The results lend support to the 
hypothesis proposed by Khorrami et al 14 that vor- 
tex shedding from the trailing edge of the slat 
is responsible for a loud, high-frequency noise ob- 
served in a corresponding set of experiments. A 2D 
FW-H solver 19 was used to evaluate the noise using 
a variety of different integration surfaces. The use 
of an on-body integration surface produced suspi- 
ciously large directivity lobes directly upstream and 
towards the microphone array. The use of off-body 
integration surfaces that include the slat-cove region 
modified the directivity pattern considerably. Alter- 
native off-body integration surfaces that include the 
slat-cove region produced similar results. Compari- 
son of mean-square pressure results with data from 
individual microphones in the acoustic array showed 
qualitatively similar trends. 

A number of 3D acoustic calculations were per- 
formed by replicating the 2D RANS data in the 
spanwise direction for varying spanwise extents. As 
expected, the noise generally increased with increas- 
ing spanwise extent. All of the important directivity 
lobes observed in the 2D acoustic results were ob- 
served at the smallest spanwise extent tested. A sim- 
ple scaling of the 3D directivity pattern compared 
favorably with the 2D directivity pattern. Unfortu- 
nately, to date, the scaling factor can only be ob- 
tained empirically, after both the 2D and 3D cases 
have been computed. A method for modeling the 
scaling factor would greatly enhance the predictive 
usefulness of the 2D solution. 
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